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Abstract 

We consider a system interacting with a chaotic thermodynamic bath. We derive an 
explicit and exact Kraus operator sum representation (OSR) for the open system 
reduced density. The OSR preserves the Hermiticity, complete positivity and norm. 
We show that it is useful as a numerical tool by testing it against exact results for a 
qubit interacting with an isolated flawed quantum computer. We also discuss some 
interesting qualitative aspects of the OSR. 
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1 Introduction 



Efforts to design scalable quantum computers (QCs) [1] and other new quan- 
tum technologies [2,3] are rekindling interest in the phenomenology of open 
quantum system dynamics. While effects of atom-radiation field interactions 
[4,5] have been extensively studied, the consequences of interaction with a 
complex environment are just beginning to be explored. Solid and condensed 
phase environments bequeath a diversity of controllable interactions for scal- 
able implementation of new quantum technologies, and thus understanding 
the open system dynamics induced by such complex environments is of par- 
ticular interest. Recent studies have reported a number of interesting effects 
such as suppression of decoherence with increasing bath chaos [6,7,8,9] and 
large coherent shifting [8,9]. 
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Manifestation of chaos [10,11], or at least non-integrability, is a defining and 
unavoidable feature of complex environments. Evidence suggests that con- 
densed phase dynamics is generally chaotic: Dynamics of a colloidal particle 
in water [12], and both vibrational dynamics [13] of Si and its electronic struc- 
ture [14] have been shown to be chaotic. Anharmonic corrections are known 
to be important in the study of phonons and essential for an understanding of 
heat transport [15]. Thus, the chaotic bath model provides a better represen- 
tation for the condensed or solid state environment than standard uncoupled 
oscillator bath models. 



Chaotic bath models are substantially more difficult to simulate than the 
radiation field. At high temperatures random matrix theories [16], or classical 
trajectory simulations such as the Wigner method [17], can be used. Redfield 
theory [18,19] and its generalizations [20,21] can also be employed. But at 
the very low temperatures relevant to quantum computing technologies the 
options are more limited. There is a recently developed approximate master 
equation [8,9,22,23,24] which appears to be accurate in this low temperature 
regime. However, as with any approximate theory its predictions need to be 
otherwise confirmed. 



In this manuscript we derive a simple formula for the reduced density which 
captures all the important features of open system dynamics induced by 
chaotic environments in the low temperature regime. The equation is of the 
Kraus Operator Sum Representation (OSR) [25] form which automatically sat- 
isfies all the required conservation laws for the reduced density, i.e. Hermitic- 
ity, complete positivity and norm conservation. Our Kraus OSR[25] provides 
a formally exact equation for the reduced density in terms of explicit represen- 
tations of the Kraus operators. Aside from chaos, the only required property 
of the bath Hamiltonian is that it be of large thermodynamic dimension. 



As a first test, we compare the predictions of the chaotic Kraus OSR[25] to 
exact calculations for a model of a single qubit interacting with an isolated 
QC with static internal flaws. While our derivation assumes a bath of ther- 
modynamic dimensions, our numerical test shows that very accurate results 
can be obtained even for a chaotic bath of just ten qubits. 



Organization of this manuscript is as follows. A derivation of the Kraus OSR[25],in 
three steps, is given in Section 2. In section 3 we discuss our chaotic spin-bath 
model. In Section 4 we test the OSR against exact numerical results. In section 
5 we discuss our results. 
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2 Derivation of chaotic Kraus OSR 



Consider the following total Hamiltonian for the bipartite closed system 

H = H S + SB + H B (1) 

where Hg is the system Hamiltonian and H B is the bath Hamiltonian. SB is 
the interaction Hamiltonian where S and B are the system and bath coupling 
operators, respectively. 

Our derivation consists of three parts. First we review the construction of the 
OSR assuming a complete bath eigenbasis (i.e H B \j) — an d J2j \ = 

Ib)- Next we show that the off-diagonal (j 7^ k) matrix elements of bath 
coupling operator 

B hk = {j\B\k) (2) 

vanish for a chaotic bath in the thermodynamic limit. Finally, we use this 
fact to show that the resulting Kraus OSR takes a particularly simple and 
potentially useful form. 

2. 1 General Kraus OSR for complete bath eigenbasis 

We consider an uncorrelated initial state for the whole system 

p(0) =p s (0)®p B (0) (3) 

where ps(0) is an arbitrary initial state of the system and ps(0) is an initial 
bath state of canonical form, i.e. 

p-Ej/kT 

Pb(0) = E— g— (4) 

The subsequent states will be given by the exact solutions 

p(t) = U(t)p(0)UHt) (5) 

of the Liouville-von Neumann equation 

dp(t)/dt = -i[H,p(t)}. (6) 
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Here U(t) = exp {— (i/h)Ht}. The exact reduced density of the system is then 
formally obtained by tracing over the bath degrees of freedom, i.e, 



p s (t) = Tr B {U(t)p(0)&(t)}. (7) 



Since the bath eigenbasis is complete, and the initial canonical bath density 
is already diagonal in eigenstates of Hb, performing the partial trace in the 
same basis yields 

Ps(t) = J2(j\U(t) (p s (0) ® ^\k)(k\) fr(t)\j) (8) 
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which in turn can be written in more compact Kraus OSR form 

Mt) = Y,^At)psmUt) (9) 



where the Kraus operators are defined by 

ic^{t) = ^r k {j\u{t)\k). (io) 

Here p k = exp {— f3E k }/Q are the initial populations of the bath eigenstates. 
Note that £j,k = Is- 

This form of the Kraus OSR is of formal interest only. Although it is exact 
and so satisfies all the required conservation laws (i.e. the Hermiticity, com- 
plete positivity, and norm) for the subsystem density it is not practical for 
computational purposes. First, the explicit forms of the Kraus operators (10) 
are not known in general, because the Kraus OSR (10) is not unique. And 
secondly, the double sum in (9) would make use of the OSR impractical even 
if simple forms for (10) were known. We will now show that in the limit of 
a chaotic bath of thermodynamic dimension, explicit forms for (10) can be 
supplied and that the double sum reduces to a unique single sum. 



2.2 Chaotic bath of thermodynamic dimension 



Using the bath eigenvalues and eigenstates, the total Hamiltonian (1) can 
alternatively be written as 

H = Hs + SY,BiAMk\+Y, E i\M\ ( n ) 

3,k j 
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where Bj^ are the bath coupling operator matrix elements defined above in 
Eq. (2), and Ej and \j) are again the bath eigenenergies and eigenstates. 

Decoherence and open system dynamical effects are mediated by the matrix 
elements Bj jk and by operators S and \j)(k\. We will focus on the properties 
of Bjj,. Now we assume that Hb is chaotic (both quantally and classically) 
and that the Wigner function [26] of B has a classical limit. Coordinate and 
momentum bath coupling operators, for example, have Wigner functions with 
classical limits. Let N denote the number of bath degrees of freedom. We will 
eventually send N — > oo to attain the thermodynamic limit. 

Quantum chaos has a number of important consequences for matrix elements 
and spectra[10,ll]. We will now show that off-diagonal matrix elements vanish 
for thermodynamic baths. We begin with an exact formula for the off-diagonal 
matrix elements of an arbitrary quantum system: 




j dy e 2^(x-x )Jy w .( x + hy/2)W k (x - hy/2). (12) 

This is obtained by combining the exact equations (50) and (119) of [27]. Here 
£>(x) is the Wigner function [26] of operator B and 

W,-(x) = Cj J dw e 2 ™ pv (q - hv/2\j)(j\q+hv/2) 

is a scaled Wigner function for eigenstate j. It is well known [28] that Hj(x) 
has a well defined classical limit if the classical Hamiltonian _f/#(x) is chaotic. 
[The Wigner function of Hb may differ from the classical Hamiltonian by small 
corrections which vanish with Planck's constant, but we ignore this distinction 
in our discussion.] We employ a notation where x = (p, q) , where p are the 
momenta and q the coordinates. J is the 2N x 2N dimensional symplectic 
matrix. Here Cj = j rfx 5(1 — Hb(^)/ Ej) and Ij^ = (CjCk) 1 ^ has the 
dimensions of action. 

Now, note that / rfx 5(1 — H b(x) / Ej) oc Eje Sj ^ kn where Sj is the microcanoni- 
cal entropy and ks is the Boltzmann constant. Now since entropy is extensive, 
Sj scales linearly with N for sufficiently large N. Consequently, by l'Hopital's 
rule, it follows that Cj scales exponentially with N. Finally, it follows that Ij^ 
is constant, or at most grows weakly with iV since the energies Ej also grow 
with N. 

The Wigner functions Wj(x) and -B(x) and the associated integrals in (12) 
have well defined classical [28] and therefore thermodynamic limits. In the clas- 
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sical limit, for example, these terms reduce to E k 5(Ej — E k )B 2 which increases 
slowly with N. (Here B 2 is the classical canonical average of the squared clas- 
sical limit of the Wigner function for B. Note that B oc y/~N and so B 2 would 
normally scale linearly with N. E k 5(Ej — E k ) is only weakly dependent on 
N. Note also that 5(Ej — E k ) should be understood as a distribution with a 
small but nonzero width.) Thus, by l'Hopital's rule, the magnitude of \Bj^\ 2 
in the thermodynamic limit is determined by the factor (-r^) N - Consequently, 

if h < Ij jk , which will always be true for sufficiently large N, then \Bj^\ 2 — > 
exponentially fast as N — > oo. Note that the requirement is h < Ij tk , not the 
stronger semiclassical h « Ij ik . 

Hence, the off-diagonal matrix elements of B vanish for a chaotic bath of 
thermodynamic dimension. They also decay exponentially fast with increasing 
N, and so they may be neglected even for quite small baths. 

In the semi-classical limit, the scaling, \Bj tk \ 2 oc h N has long been known in 
the quantum chaos literature [27,29], but its consequences in the thermody- 
namic limit, and for the Kraus OSR in particular, have not been appreciated. 
As we have shown, the semi-classical limit h — > is not a necessary condi- 
tion for the vanishing of the off-diagonals, \Bj :k \ 2 — > 0, for a chaotic bath of 
thermodynamic dimension. 

It is also worth noting that similar formulas arise in other contexts from com- 
pletely different arguments [30]. 



2.3 Chaotic Kraus decomposition 



Now, we will show that the vanishing of the off-diagonal matrix elements of 
the bath coupling operator has important consequences for the Kraus OSR 
(i.e. Eq. (9)). To begin with, Hamiltonian (11) reduces the simpler form 

H = H S + J2(SB k , k + E k )\k)(k\. (13) 
k 



For integer powers of I — 0, .., oo then, it is true that 

(j\H l \k) = (H s + SB ktk + E k ) l 6 hk (14) 

for all j and k. Using Eq. (14) for 

e -i* = £MW (15) 
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one can readily show that the Kraus operators (10) take the explicit form 

= y/pie-W* +aB ** +E #6 jlk . (16) 



Finally, substituting this back into (9) gives the desired form of the chaotic 
Kraus decomposition 

p s {t) = Y, e ——e-^ Hs+SB ^ t p s {Q)e^ Hs+SB ^ t . (17) 

k Q 

Equation (17) is an exact equation for the reduced density of an arbitrary 
quantum system interacting with a chaotic bath of thermodynamic dimension. 
While the arguments of the previous section were cast in terms of a phase 
space representation, similar arguments could be devised for spin baths using 
generalized coherent state representations [23]. In the next section we actually 
employ a spin bath to test the formula (17) numerically. We will show that 
(17) is very accurate even for a bath of just ten qubits. 



3 Test model 

In our recent investigation [8] we showed that internal decoherence, dissipation 
and coherent shifting, caused by static one and two-body internal imperfec- 
tions, are potential error sources for flawed QCs[31]. We also showed that a 
single detector qubit can be configured to probe internal bath dynamics [9] 
and thus one can obtain valuable information about the bath self-interactions. 
Here we extend our previous exact numerical results to a longer time limit and 
test the chaotic Kraus decomposition against these exact results. In what fol- 
lows we briefly review our flawed QC model and our exact numerical approach. 
Details of this study are given in [9] . 

3.1 Exact numerical approach 

Our chaotic spin-bath model represents a flawed QC core[31]. The total Hamil- 
tonian takes the form 

H = -\B&?>+*V> (18) 

1 N N-1 N 

- ± E Of ^ + Blaf) + E E J?*S>*y, 

Z i=l i=l j=i+l 
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where the first term represents a central qubit (i.e. the detector qubit). The 
second term is the system-bath coupling operator, and the last two terms 
constitute the chaotic bath Hamiltonian which represents the rest of the flawed 
QC. 

Our numerical simulations are based on the experimentally relevant parameter 
regime for a Josephson charge-qubit QC proposal [32]. We set the detector 
qubit energy to Bq = 1 e in units of e = 200 mK. We assume that all other 
qubits other than the detector qubit have one-body flaws, i.e. Bf e [Bq — 
5/2, Bl + 5/2] and similarly Bf £ [B% - 5/2, B* + 5/2] where the detuning 
is set to 5 = 0.4 e. Two-body flaws, i.e. residual system-bath couplings, and 
chaos generating intra-bath interactions are modeled via A, G [—A, A] and 
J*'- 7 G [—J xi Jx], respectively. We kept the bath temperature constant at 
kT = 0.25 e in all our simulations. Keeping the value of the system-bath 
coupling constant at A = 0.05 e we considered three different values of intra- 
bath couplings; J x = 0.50, 1.00, 2.00 in units of e. By J x = 2e the bath is fully 
chaotic. 

Exploiting the experimentally relevant low temperature regime [33] for the 
charge-qubit QC proposal [32] we calculated the exact reduced density ps(t) 
via the formula 

20 

Ps(t) = T,Pn Tr B[\* n (t)){* n (t)\], (19) 
n=l 

where the populations of bath are defined by 

e -E n /kT 

Pn = V 20 P-E m /kT ( 20 ) 



and \^ n (t)) are the solutions to the Schrodinger equation for all initial states 
| x I / n(0)) = | V(0)) <8> \n). The initial system state was chosen as a superposition 
state |^(0)) = (|0) + \ l))/y/2. 20 bath eigenstates sufficed for all our calcula- 
tions at the given low temperature. We performed numerical integrations by 
an eighth order variable stepsize Runge-Kutta method [34] and exact diago- 
nalization of the bath Hamiltonian is achieved by a Lanczos algorithm[35] for 
N =10 bath qubits. 

3.2 Chaotic Kraus decomposition approach 

The explicit form of Kraus operators for the Hamiltonian of our model is 
iC n (t) = V^e^^-^^)*. (21) 
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The chaotic Kraus decomposition for the given temperature kT and initial 
density of the form p(0) = \i/j(0))(i/j(0)\ then has the form 



where 



20 



Ps(t) = Y,Pn 

n=l 



\cUt)\ 2 cUt)[c n (t)V 

c%(t)[cUt)r K{t)\ 2 



(22) 



<?(*) : 

c n o(t) 



V2 
2 

V2 
2 



\B Z — 2_B n n ) 
cos a n t + i ; '■ — sin a n t 

On 

. (B z + 2B nn ) . 
cos a n t — % sin a n t 

On 



(23) 
(24) 



and b n = ^Bl + 4B~l n 



(25) 



We calculated the diagonal matrix elements of the bath coupling operator by 
using the exact eigenstates of the bath Hamiltonian. 



3.3 Test for chaos 



In order to confirm an integrable to chaotic transition in the bath Hamilto- 
nian Hb as a function of J x , we analyzed the nearest neighbor level spacing 
distribution[10] P(s) for the 200 unfolded lowest eigenvalues. The unfolded 
spectrum Ei is obtained from the actual eigenspectrum Ei via smoothed stair- 
case functions N(E), i.e. Ei = N(Ei). Smoothed staircase functions are poly- 
nomial best fits to the actual staircase functions N(E) = X^=| Q(E — Ej) 
where Q(x) is a Heaviside step function. The onset of chaos is indicated by 
the functional form of P(s) changing from the Poisson form exp (— s) charac- 
teristic of integrable systems to the Wigner-Dyson form (7r/2)sexp (— 7rs 2 /4) 
associated with chaos. Our results are given in Figure 1. The onset of chaos 
is observed above J x = 0.15 e. For stronger couplings the eigenstatistics are 
consistent with Wigner-Dyson distribution. 

While this level spacing statistic is considered to be a universal indicator of 
quantum chaos, it does not provide information on the degree of chaoticity. 
We wish to know how the accuracy of the Kraus decomposition changes with 
increasing degree of chaos. Therefore, we alternatively refer to the Loschmidt 
echo M(t)[36], which is also widely believed to be an efficient indicator of 
quantum chaos [3 7]. 
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Fig. 1. Level spacing statistics for chaotic bath Hamiltonian 
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Fig. 2. (a) Short and (b) long time Loschmidt echo dynamics of bath Hamiltonian 
for J x = 0.50 (black solid line), J x = 1.00 (red dotted line), and J x = 2.00 (green 
dashed line). 

We calculate M(t) for the bath Hamiltonian with the following formula 

M(t) = \^ \exp{tH Q t/h}exp{-t(H + V)t/h}\^ )\ 2 (26) 

where \ipo) is an initial state chosen as the ground state of Ho. H Q is the 
integrable bath Hamiltonian (i.e. Hb for J x = 0.00) and V is the chaos gener- 
ating perturbation Hamiltonian for J x = 0.50, 1.00, 2.00. A summary of M(t) 
calculations is presented in Fig. (2). It is clear from Fig. (2) that increasing 
the magnitude of intra-bath coupling J x results in faster exponential decay of 
M(t), which can be interpreted as increasing the degree of chaoticity. 
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3.4 Quantitative measures of environmental effects 

We use purity, V(t) = Tr s [ps(t)} 2 , to determine the extent of non-unitary 
evolution in the system (i.e. decoherence and dissipation). The ideal value of 
purity is unity for pure initial system states. The fidelity is defined by J-{t) = 
Tr s[f>s{t)p l s eal {t)} where p^ eal (t) is ideal system evolution in the absence of 
system-bath interactions. The ideal value of fidelity is also 1. While the purity 
is insensitive to the unitary effects of system-bath interactions, the fidelity is 
sensitive to both unitary and non-unitary effects. Therefore, large deviations 
between the magnitudes of purity and fidelity are a good indicator of the 
presence of unitary shifts. 



4 Results 

We tested the accuracy of the chaotic Kraus OSR for several values of intra- 
bath couplings in the chaotic bath regime. Here we compare the prediction of 
the chaotic Kraus decomposition with exact numerical results by calculating 
the purity V(t) and fidelity F{t) of the reduced density. 

In Figure (3) we plot V(t) for three different values of intra-bath coupling J x . 
(Time is in units of H/e.) The dashed lines are the Kraus predictions while the 
solid curve represents the exact dynamics. Even for J x = 0.5 e, shown in part 
(a), the decoherence predicted by the Kraus decomposition is of the correct 
order of magnitude. Clearly, however it is not quantitatively accurate. There 
is also a faster time scale to the exact dynamics which is not captured at all by 
the Kraus decomposition. In part (b), for J x — 1 e and hence stronger chaos, we 
see much better agreement. In part (c), for J x = 2 e, the Kraus decomposition 
is in very good agreement with the exact results. This is remarkable given 
that the bath contains only ten qubits. Moreover, it suggests that the Kraus 
decomposition could actually be employed as a useful computational method 
for low temperature systems. 

In Figure (4) we plot T{t) for both exact and Kraus decomposition results 
for the same three different values of intra-bath coupling J x . Notice that the 
magnitude of the errors in J-{t) is much larger than those in the purity V(t). 
The purity measures only non-unitary errors. J- if) is also sensitive to unitary 
errors. The large deviation of Vit) from unity therefore indicates the presence 
of a large coherent shift. The origin of this shift has been explained in detail 
elsewhere [9]. The agreement between the Kraus and exact results is quite good 
even for J x = 0.5 e (shown in part (a)). For J x — 1 e (part (b)) and J x = 2 e 
(part (c)), the agreement is excellent. 



11 




(c) 




Fig. 3. Exact numerical (solid lines) and Kraus (dotted lines) results for purity V(t): 
for J x = 0.50 in (a), J x = 1.00 in (b), J x = 2.00 in (c). 
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Fig. 4. Exact numerical (solid lines) and Kraus (dotted lines) results for fidelity 
F(t): for J x = 0.50 in (a), J x = 1.00 in (b), J x = 2.00 in (c). 
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5 Discussion 



The Kraus decomposition captures the correct quantitative behavior for both 
unitary and non-unitary types of contributions from chaotic environmental 
dynamics. The accuracy of the Kraus decomposition even for our small bath 
model is quite impressive and thus it it may prove a useful numerical simu- 
lation method. The suggestive form of the decomposition (17) is the result of 
primary interest. For example, it suggests that all non-unitary effects vanish 
when all diagonal coupling matrix elements are identical, i.e. B k k = B. In this 
case the dynamics reduces to 

p s (t) = e -^s+SB)t h ^yi { H s+ SB)t (27) 



and the only effect of the bath is the coherent shift H s — > H s + SB. 

Moreover, it is also clear that the greater the variance in the coupling ma- 
trix elements, the more decoherence will be observed. This prediction is in 
agreement with that of an approximate master equation[8,9,22,23,24] which 
has been shown to be quantitatively accurate for a number of systems[8,9]. 

It also predicts that a chaotic bath of thermodynamic dimension with uncor- 
rected initial state cannot cause decoherence or dissipation at absolute zero 
temperature. Such a bath can still cause coherent shifting, however. It was 
previously known that strongly chaotic finite baths cause little decoherence at 
low temperature[6,7]. The fact that no decoherence is possible at zero temper- 
ature for large chaotic baths is in strong contrast to the behavior of integrable 
baths. In the spin-boson model strong decoherence is possible even when the 
bath has zero temperature. 



6 Summary 



We derived an exact and explicit Kraus decomposition for systems interacting 
with chaotic baths of thermodynamic dimension. We showed that the decom- 
position is very accurate even for a small finite qubit bath, and it predicts all 
important features of the open system dynamics induced by chaotic environ- 
ments. It could therefore prove useful as a numerical simulation method. 
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